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Abstract 



^D I Beginning with Bekenstein, many authors have considered a uniformly spaced 

^— V ' discrete quantum spectrum for black hole horizon area. It is also believed that 

CN ■ the huge degeneracy of these area levels corresponds to the notion of black 

^^ . hole entropy. Starting from these two assumptions we here infer the algebra 

of a Schwarzschild black hole's observables. This algebra then serves as moti- 
vation for introducing in the system's Hamiltonian an interaction term. The 
interaction contains the horizon area operator, which is a number operator, 
and its canonical conjugate, the phase operator. The Hawking radiation from 
^ . a Schwarzschild black hole is seen to be a consequence of an area-phase inter- 

^«0 I action. Using this interaction we have reproduced the semi-classical result for 

the Hawking radiation power. Furthermore, we show that the initial state of 
^O ' the black hole determines the nature of its development. Thus, a state which 

^ ■ is an area eigenstate describes a static eternal black hole, but a coherent state 

Q^ . describes a radiating black hole. Hence, it is the observer's initial knowledge 

CJ ! or uncertainty about the horizon area which determines the evolution. 
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I. INTRODUCTION 

The first quarter of tliis century yielded two revolutionary theories, quantum mechanics 
and general relativity; they have changed drastically the way we think about the physical 
world. Since then many efforts have been made to unify these theories, and today we are 
presented with a variety of generic "quantum gravity" theories. Let us mention canonical 
quantum gravity [l|], perturbative quantum gravity 0, superstring and D-brane theory 0, 
loop (nonperturbative) quantum gravity Q] and others. All these are intricate theories: even 
comparison of their predictions is a difficult task. In view of the dearth of concrete experi- 
mental results bearing on the quantum aspects of gravity, one wonders whether it would not 
be more to the point to construct a theory which emphasizes some particular gravitational 
system where circumstances allow a guess as to the general nature of quantum effects. Black 
holes suggest themselves as a good place to start: they involve strong gravitational fields, 
have properties reminiscent of localized particles, and results from quantum field theory 
[principally Hawking's radiance (HR)] suggest that quantum effects are central. Thus an 
approach to black hole quantization has emerged 0-§l which starts off from hints garnered 
in classical black hole physics, and proceeds to quantum ideas. 



This approach hails back to the realization by Christodoulou |^, Hawking [|I0[ and Pen- 
rose and Floyd [^ that transformations of a black hole (BH) generically have an irreversible 
character. One of its conclusions is that the horizon area tends to grow, a rule which served 
as motivation for Bekenstein [0 to associate entropy with horizon area A according to (we 
insert the modern value of the coefficient) 

S = A/{4Cl) (1.1) 

(£p and A4p will denote the Planck length and mass, respectively). Christodoulou also 
showed that there are some special BH reversible processes which leave BH area unchanged. 
This has led Bekenstein [Q to remark that the horizon area of a non-extremal BH behaves 
as a classical adiabatic invariant (for support of this conjecture see also Bekenstein |TB 



and Mayo [T^). According to Ehrenfest's principle |jT5[| any classical adiabatic invariant 



corresponds to a quantum entity with discrete spectrum. Thus Bekenstein conjectured that 
the spectrum of the horizon area of a (non-extremal) BH should be quantized. Today the 
idea of a discrete eigenvalue spectrum for the horizon area is also supported by the work of 
Ashtekhar and others Q]. 

In support of a conjectured uniform spacing of the area eigenvalues, Bekenstein noted P,|^ 
that the assimilation of a neutral test particle into a BH increases the horizon area by at 
least vCj:,, where t; is a dimensionless constant of order unity. Also Hod |l^ have shown 



independently that the same kind of lower bound on the horizon area increase applies to 
the BH assimilation of a charged particle. The fact that the lower bound turns out to be 
independent of the BH parameters (mass, charge and angular momentum) is in harmony 
with a uniformly spaced area spectrum: 

an = vCl{n + r]); < r] < 1; n = 0,1,2,... (1.2) 

Here a„ is the area eigenvalue. In the present work we shall set r^ = for simplicity. 

The spectrum in Eq. ( |1.2| ) together with the classical relation A = IQttM^ (we use 
gravitational units in which G = c = 1) implies that for n ^ 1 the energy (mass) levels M„ 



of an isolated Schwarzschild BH have the form M„ = moy/n; here mo = Jvh/16n. This 
kind of spectrum is compatible with Bohr's correspondence principle. That is, for large n 
values the n -^ n — 1 transition frequency 

^^ = h = 32^ ^'-^^ 

agrees with the classical BH oscillation frequencies which scale as M~^. 

The energy Uujq emitted in the transition n -^ n — 1 is reminiscent of the Hawking temper- 
ature. Hawking |r^ has shown that a steady emission of particles to infinity originating in 



pair creation in the gravitational field follows when a BH is formed by gravitational collapse. 
He associated with the BH a temperature Tbh = h/8nM [compatible with the entropy given 
in Eq. (|1.1| )] which characterizes the thermal spectrum of the emitted particles. Note that 
Tbh and huJo are of the same order. Hence, the idea that the temperature represent the 
average energy of the emitted particles is compatible with a mass spectrum M ~ ^/n. 

Following Bekenstein's conjecture about the mass spectrum, a number of authors (see 



full references in [iH] and M) have obtained the selfsame spectrum from a wide variety 



of arguments. For example, Vaz and Witten |]T9[ examined the Wheeler-DeWitt equation 



for a static, eternal Schwarzschild BH in Kuchaf-Brown variables, and obtained the mass 
spectrum M„ oc ^/n. Louko and Makela [^] have proposed a Hamiltonian operator for 
the Schwarzschild BH, and obtained the same spectrum. Kastrup |^I| explained how such 



a spectrum may be understood in the framework of canonical quantization of a purely 
gravitational spherically symmetric system. Carlip and Teitelboim IQ used the standard 
(Euclidean) action principle for the gravitational field to show that the horizon area and the 
opening angle at the horizon are canonical conjugates. The opening angle is bounded (in 
the Euclidean sector). This is compatible with the idea that the area operator is represented 



by a number operator (see Eq. ([T7^ ), since its canonical conjugate is a phase operator [^ 
This similarity is further evidence in favor of an equally spaced area spectrum. 

The energy levels M„ are expected to be degenerate. Denoting the degeneracy by g{n) 
and identifying In g{n) with the entropy given in Eq. ( [1 . 1| ) , Bekenstein and Mukhanov ^ 
found that i; = 4 In A; [t; is the dimensionless constant in Eq. (|1.2|) ], or equivalently g{n) = /c" 
where /c = 2, 3, 4, ■ ■ • They chose k = 2 for simplicity; recently Hod [^ has given evidence 
in favor of the choice k = 3. 

By contrast to the approach just discussed, quantum loop gravity determines a non- 
uniformly spaced area spectrum. Ashtekhar and Krasnov claim [Q that the area spacing 
approaches zero exponentially for large quantum numbers. Hence, the area spectrum they 
propose is directly compatible with the semi-classical prediction by Hawking |T^ that the 



emitted spectrum is smooth. The present paper shows that other features of Hawking radi- 
ation emerge naturally from the uniformly spaced area spectrum hypothesis. For example, 
we show that the rate of change of the mass of a BH (not near the Planck scale) agrees with 
that of the semi-classical calculation by Page |]2S| up to the second order in time. Thus, 



although the uniformly spaced area spectrum is inimical to the smooth Hawking spectrum, 
it is not at variance with equally important features of a semi-classical BH. 

It is clear that the Hawking temperature is a quantum effect {Tbh is proportional to 
Planck's constant h), and a deep understanding of it requires a quantum theory of BHs. In 
this paper a quantum model for Schwarzschild BHs is proposed. We start from two axioms: 



(a) that the horizon area spectrum is equally spaced, and (b) that the degeneracy of level n 
is A;". Surprisingly, these two axioms enable us to obtain the algebra of the Schwarzschild BH 
observables, and even the form of the Hamiltonian. The algebra obtained is the algebra of 
angular momentum with the restriction that the absolute value of the "angular momentum" , 
here called the hyperspin, can have only the values (fc" — l)/2. This last limitation motivates 
us to introduce an interaction term which allows transitions only between allowed values of 
the hyperspin. 

The z-projection of the hyperspin represents the internal degree of freedom of the BH 
[this degree of freedom corresponds to the entropy given in Eq. ( p..!! )]. According to the 
"no hair" theorems ||2^, an observer located at large distance from the BH is not able to 



distinguish between different states of the internal degrees of freedom. Hence, we find it 
necessary to introduce a horizon projection operator which filters out the hyperspin degree 
of freedom, but leaves that corresponding to the mass of the BH. This corresponds to the 
conception that most of the information about the BH interior is hidden (by the horizon) 
from an exterior observer. 

The Hilbert space which describes the states that the observer at infinity can measure is 
spanned by area (mass) eigenstates. We show in this paper that the initial state determines 
the nature of the dynamics of the BH. For example, if the initial state is chosen to be an 
area eigenstate, it describes a static eternal BH. On the contrary, if the initial state is chosen 
to be a coherent state,[] it reproduces the semi-classical result l2^J3l ^^^ ^^^ Hawking power. 
That is, it is our initial knowledge (uncertainty) about the horizon area which determines 
the nature of the dynamics. A BH formed by gravitational collapse generally has uncertain 
horizon area, and thus it looses mass, whereas an eternal BH, with a definite mass and area, 
is static. 

This paper is organized as follows. In section [11 A] we construct the hyperspin algebra 



and characterize the observable Hilbert space. In section [II we define the mass operator. 



and then find the general form of the Hamiltonian of a Schwarzschild BH. We also review 



the quantum phase problem in general, and then define the phase of a BH. In section ^ 
we discuss the initial states to be considered, and in section we use them to compare 
semiclassical properties of the HR with our model's predictions. Finally, in section |V11| , we 
present our summary and conclusions. 

II. BLACK HOLE ALGEBRA 

The classical Schwarzschild BH is characterized by only one parameter, its mass or 
equivalently its horizon area. Since the horizon area is an adiabatic invariant and is suggested 
to be quantized (as discussed in the Introduction), we shall focus on it. Hence, the first step 
should be to investigate the horizon area operator. In this work we shall not construct the 
area operator from more fundamental quantities |Q, even though this might illuminate our 
perception of the internal structure of the BH. Instead, we shall use the properties of the 



^Coherent states are the ones that most closely imitate the classical results for a simple harmonic 
oscillator. 



area operator which follow from general arguments (see the Introduction). We now turn to 
construct the algebra. 

A. The Algebra 

Let A represent the horizon area operator (boldface is used everywhere to denote oper- 
ators). Motivated by the arguments discussed in the Introduction we shall start from the 
following two axioms: 

(1) The eigenvalues of A are a„ = aon where n = 0,1, 2..., ao = 41n/c£p is a positive univer- 
sal constant, and k is an integer greater then one. 

(2) There are exactly fc" independent eigenstates of A with eigenvalue a„, i.e. the degeneracy 
of a state with area a„ is k^. 

From the second axiom it follows that there is at least one more operator which commutes 
with A. Let G represent the set of all independent operators which commute with A. 
Typically, 

[A,G] = 0, (2.1) 

and the set {A, G} of "observables" is maximal in the sense that if S is another observable 
which commutes with both A and G, then S = /(A, G) where / is a function of A and G. 
We shall call G a "secret" observable since it describes the internal degrees of freedom of 
the BH. 

From the second axiom or from Eq. ( p.l| ) it follows that we can simultaneously diagonalize 
A and G. Hence we shall denote the states as \n,m) such that 

1^^171,771) = n\n,m) ra = 0,1,2... 

G\7i,7n) = 7n\7i,7n) tti = 1,2, ...,k"' (2-2) 

where N = aQ^A is a dimensionless number operator. The eigenvalues m of G are limited 
by n since we want a degeneracy of fc". We have assumed that tu is an integer; if not, instead 
of G we shall choose the "secret" operator to be /(G), where f{m) = 1, 2, ..., fc" for all the 
different fc" values of 7n. 

Now, since G has a discrete spectrum m = 1,2, ..., k"-, we shall introduce the annihilation 
operator g and the creation operator g^ corresponding to G: 

[G,g] = -g [G,gt]=gt. (2.3) 

As we shall see shortly, under certain conditions g and g^ have to satisfy a certain commu- 
tation relation. In order to single out an algebra, we have to demand that [g, g^] be linear 
in G. In order to derive this commutation rule we shall first guess some of the properties 
of g and g^ from physical arguments. First we do not want g and g^ to change n when 
they act on states with a definite area eigenvalue n. This means that transitions among the 
internal degrees of freedom {m — > m') do not have to change the area of the BH. Thus, the 
area operator N commutes with g and g^ 

Next, from Eq. ( |2.3| ) we have by analogy with the harmonic oscillator 

g^\n,m) (X \n,7n + 1). (2-4) 



Therefore g^|n, A;") = and thus also gg^|n, fc") = 0. It is clear from Eq. ( |2.3| ) that [G, gg^] = 
and also that [N, gg"''] = 0. But the set (N, G) is maximal and thus, since gg''" is an 
observable, gg^ = /ii(N, G) where hi is a function of A and G. Hence, from the fact that 
gg^|n. A;") = we find that /ii(n, /c") = 0. Similarly we find that 

gtg|n,l) = 0, (2.5) 

where g^g = /i2(N, G). It follows from Eq. ( p.5|) that h2{n^ 1) = 0. Collecting all of these 
facts we find that gg''" and g^g can be written uniquely as 

ggt = m(N)G(A;N - G) 

gtg = M(N)(A;N_G + l)(G-l), (2.6) 

where m(N) is some function of N. In Appendix A we show that m(N) = 1. Thus, we settle 
on the commutation relation 

[g,gt] = A;N-2G + l (2.7) 

which is linear in G. Eq. (|2.71 ) and Eq. ( p.3|) seem to imply that the set of operators 
{G,g,g^} plus the identity element span a 4-dimensional Lie algebra which constitutes a 
central extension of the rotation group. The representations are supposed to be such that 
the maximal value of the third component is just A;". The type of the central extension is 
exhibited in Eq. ( |2.7| ) where the term (fc" + l)/2 is the central charge (a multiple of the 
identity). However, as we shall see later, we can eliminate the central charge from Eq. ( p.7|) 
by a redefinition of the third generator G ^ G — {k^ + l)/2. Furthermore, we shall make it 
clear under what circumstances the commutation rule ( p.7|) is consistent with the required 
spectrum of N and G. 

1. The Matrix Elements 

We turn now to calculate the matrix elements of the operators in the algebra. We chose 
the orthonormal set of states {|ra, m)}^"lj' 2 '■fcn as a basis. Thus we have 

(nW|N|n, m) = n5n,n'5m,m' (2.8) 

(nW|G|n, m) = rn5n,n'5m,m'- (2.9) 

Now, from (|2.4|) it follows that g"'"|n, m) = C^^\n,m + 1) where 



l^nmP — {nm\g,'g^\n, m) = {n'm\G{k^ — G)\n, m) = m{k^ — m). (2-10) 

In a similar fashion g|n, m) = C^^\n, m — 1) were 

l^nmP — {nm\g^ g\n, m) = {nmKk^ — G + 1)(G — l)\n, m) 

= (A;"-m + l)(m-l) (2.11) 

and thus we finally obtain for the matrix elements of g and g^: 

{n'm'\g^\n,m) = {m{k'^ - m))^'''^ 5n,n'Sm',m+i 

{nm'\g\n,m) = {{k'' - m + l){m - l))^^"^ 5n,n'Sm',m-i, (2.12) 

where for simplicity we have chosen C^^ and C~^ to be a positive real numbers. These 
matrix elements are reminiscent of the matrix elements of J_ = J^ — iJy and J_ = J^ + iJy 
in the algebra of angular momentum. As we shall see shortly, this is not accidental. 
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2. The Hyperspin of a Black Hole 
Using Eq. ( |2.(j| ) we can find a relation between N, g and g''" 

fc''' = [g, g? + 2{g, gn + 1 (2.13) 

where the symbol {, } indicates an anticommutator. It is not clear a priori that building N 
from g and g^ in this way is consistent with the required spectrum of N. In order to clarify 
this point let us first define three observables: 



Si 


= ^is + s^) 




S2 


= ^is-s^) 




Ss 


.G-i.-- 


1 

'2 



(2.14) 

Using the commutation relations between g, g\ N and G we notice that the operators 
defined in ( p.l4 ) satisfy the standard commutation rules of an angular momentum, that is 



[Si,Sj] — iSijkSk, (2-15) 

where eijk is the standard antisymmetric tensor. This angular momentum algebra is not 
surprising for an algebra with three generators. As we shall see below, the requirement that 
N is a number operator complicates the algebra. 

Let us now rewrite Eq. ( p.l3| ) in terms of the angular momentum operators: 



,2N 



48^ + 1 (2.16) 



where S^ = S^ + Sg + S| = s{s + 1) (s is a half integer). It is clear from Eq. ( |2.16| ) and the 
definition of S3, that a simultaneous eigenstate of S^ and S3 , which is denoted by |s,m3), 
is also a simultaneous eigenstate of N and G. Explicitly, |s,m3) = \n,m) for 

m = s + ms + l and P = 2s + 1. (2.17) 

Two features follow from this relation. Firstly, the modulus of the angular momentum 
(the "hyperspin" of the BH) is a true observable, whereas its 2;-projection represents the 
hidden degrees of freedom of the BH. Secondly, we can see that n may be an integer, but 
it may also be logfc(2s + 1) where s = 0, |, 1, |,2, .... Thus, the Hilbert space H spanned 
by the set {\n, m)} for an integer n is a subspace of the "full" Hilbert space H{ spanned by 
the set {\s,m^)}. Since we have obtained the angular momentum operators, one can build 
the algebra in the opposite direction (starting with the three angular momentum operators) 
and ask, what is the physical reason that a BH state is restricted to a proper subspace of 
Tif? Why can it not have general values of s, but only those of the form (/c" — l)/2? 

These questions motivate us to enter into the dynamics, and define the Hamiltonian. It 
is physically clear that transitions between states with spin of the form (A;" — l)/2 and states 
which are not of the same form are forbidden. Before starting to construct the Hamiltonian, 
let us introduce the Hilbert space of an observer who is located very far a way from the BH. 



B. The Observable Hilbert Space 

An observer far away from a Schwarzschild BH is able to measure only the mass of the 
BH. "No hair" theorems |^ suggest that this is the sole parameter characterizing the BH. 
On the other hand, the states in Eq. (|2.2| ) are characterized by the two parameters n and m. 
This raises the question, how can we harmonize between our states and "no hair" theorems? 
In other words, how to include the Israel-Carter conjecture that "black holes have no hair" 
in the Hilbert space which is spanned by the degenerate area eigenstates? 

At first sight, a description of BHs by means of mixed states with entropy S = A/AC'i in- 
stead of pure states (zero entropy) seems to give the answer to the question above. Denoting 
the density matrix by p„ for each area eigenvalue n, we find that 

Pn = T;i^\n,m){n,m\ (2.18) 

'^ m=l 

and thus the entropy S = — Trp„lnp„ = A/AC^. This agrees with Hawking's principle 



of ignorance p9[: "the surface emits with equal probability all configurations of particles 
compatible with the observers limited knowledge" . However, the problem with Eq. ( |2.18| ) is 
that it suggests that the observer is able to identify \n, mi) as a different state from \n, m2) 



where mi ^ m2- Thus, Eq. ( |2.18| ) is not compatible with "no hair" theorems. But without 



a density matrix, how to attribute a BH an entropy S = A/AC'il In order to resolve this 
issue we introduce the concept of a secret observable. 

Motivated by "no hair" theorems we claim that the observer is not able to distinguish 
between different states with the same area eigenvalue. Hence there is a subspace of Ti which 
represents the observable Hilbert space. In this subspace for each area eigenvalue there is 
only one area eigenstate, i.e. the area operator is not degenerate in the observable subspace. 
Now, each area eigenstate in the observable subspace is a linear superposition of the states 
{|n,m)}m=i,2,...,fc", that is 

\n)ob = E Cl^, m) where ^ ICI' = 1- (2-19) 

At first sight Eq. (|2.19| ) implies that an observer at spatial infinity would have to attribute 
the BH a vanishing (Bekenstein) entropy. However, the observer is not able to measure 
the coefficient (i!J^, and this lack of information corresponds to the Bekenstein entropy. The 
observable Schwarzschild Hilbert space Tiob is then defined as the subspace of Ti which is 
spanned by the states defined in Eq. (|2.19|) : 

Hob = span{\n) oh] ■ (2.20) 

We shall now explain why the above definitions are relevant to the question discussed 
at the beginning of this section. The horizon splits the space in two parts. All the physical 
processes which occur in the interior part of the BH spacetime can be viewed from the 
exterior only through a change in the mass of the BH. Hence, physics inside the horizon 
should be described in terms of the states |n, m) G 7i, whereas the physics as it is viewed 
at infinity should be described in terms of the observable states \n)oh G 'Hob- That is, the 
horizon filters out the knowledge about the internal degrees of freedom. 
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We shall define the horizon projection operator: 

oo 

h=Y,\n){n\ (2.21) 

ra=0 

where from now on \n) = \n)ob- Note that the horizon projection operator is the unit operator 
in the observable subspace. The projection operator represents the horizon in which all the 
parameters describing the BH are hidden apart from its mass. That is, each state \tp) in H 
and each observable B in 7i are reduced to 

B — >Bob = hBh (2.22) 

when a physical measurement is made by an observer located at large distance from the BH. 
Eq. ( |2.22|) is a manifestation of the concept of a secret observable. It can be regarded as 
a quantum version of the "no hair" theorems because it implies that an observer at infinity 
cannot predict the internal state of the BH apart from its mass, or equivalently its area (the 
horizon projection operator erases the knowledge about G). There are, indeed, two systems 
involved. One is the "real" physical system of the BH (without the horizon filter) which is 



described by the mixed state given in Eq. ( |2.18|) . Due to the internal degrees of freedom the 
system has an entropy S = A/ACp. The other system is a small 'slice' of the real physical 
system ( p.2(]| ), and only this 'slice' is accessible to the observer at infinity. One can view it 
as if the observer is wearing a pair of glasses that filter out all but one color (that is, the 
observer state \n)ob)- We now turn to find the form of the observable Hamiltonian. 

III. THE HAMILTONIAN 

There are two properties which an observer located at large distance from the BH can 
measure: the mass of the BH (including the energy in the medium between the BH and the 
observer) and the energy flux (if there is one) coming from the BH. Hence, we shall start 
this section by defining an operator in the algebra which represents the mass of the BH. 

A. The Mass Operator 



In section P A| the set (A, G) had been chosen to be maximal in the sense that if S is 
another observable which commutes with both A and G, then S is a function of A and G. 
Now, in the observable Hilbert space, the above argument reduces to the following one: if 
S is an observable in Tiob and [A, S] =0 then S is a function of A. We shall assume that 
the mass operator M which represents the BH mass commutes with A. That is, M is a 
function of A. 

Classically, the mass of the BH can be obtained from its area by the well known relation 
A = IGvrM^. We shall assume that to a good approximation this relation also holds true for 
the operators, and hence the mass operator can be written as 



M^jA(i + ,(A)) (3.1) 



where the eigenvalues of M constitute the mass spectrum of the BH. The dimensionless 
function e(A) must approaches zero in the hmit of a massive BH in order to reproduce the 
classical relation A = IGvrM^. Now, since e(A) is a dimensionless function of A it can be 
approximated, to first order, by some power of £p/A. Hence, its contribution to Eq. ( p.l|) 
may be significant only for BHs of the order of the Planck mass. We shall be interested in 
BHs not near the Planck mass, and thus we redefine the mass operator as 



M = J— (3.2) 

V 167r ^ ^ 



B. General Form of the Hamiltonian 



We shall now try to find the form of the Hamiltonian. If we assume that it is equal to 
M, we can not understand why the hyperspin of the BH can have eigenvalues only of the 
form (A;" — l)/2. Moreover, this Hamiltonian describes a static system which clashes with 
the existence of the HR. 

Let us now list the demands from the correct Hamiltonian: 
(1) It should be a Hermitian operator. (2) For a very large mass the Hamiltonian of the 
system should be equal to the mass of the black hole. (3) It should allow transitions only 
between states with hyperspins s and s' which satisfy the condition s' = s ± k^/2 for 
some integer n. (4) It does not contain G. The last demand arises because in the observer's 
Hilbert space the area operator is not degenerate (see section [II B| ). From these requirements 



it follows that we can write the observable Hamiltonian in the form 

H = M + V (3.3) 

where V describes an "interaction term" which approaches zero as the BH mass becomes 
very large. 

In section |I| we have defined three Hilbert spaces. The first one, is the "full" Hilbert 
space Hf which is spanned by the set of hyperspin states {\s, m^)}. The second is a subspace 
Ti of the full Hilbert space which is spanned by the set of states {|n, m)} where n and m are 
given in Eq. (|2.17| ) . The third Hilbert space is a subspace Hob of TC which is spanned by the 
observable states defined in Eq. (|2.19| ). In section [II B| we have explained why, in the eyes of 
an observer at infinity, Tiob is the relevant space to work with. The question remains, why 
we did have to define Hob as a subspace of H (with area spectrum a„ = aon) and not of 
Hf (with area spectrum a^ = log;, (2s + 1))? Equivalently, why should the area eigenvalues 
be equally spaced? The answer to the questions above lies in the third requirement on 
H. Let |s) and |s') be some states in Hf with hyperspin s and s' respectively. Then, by 
implementing the third requirement on H we find that (s'|H|s) = unless s = s' ± A;"*/2 
where n* is an integer. Thus, if the initial state of the BH, |s), belongs also to H (that is 
s = {k^ — l)/2), then (s'|H|s) = unless \s') G H; hence the evolution of the system takes 
place in a subspace H of Hf. 

On the other hand, if the initial state \s) does not belong to H then (s'|H|s) = for any 
|s') eH so that H does not play any role in the evolution of the BH. In this case, the area 
eigenvalue of the initial state \s) is a^ = logj;,(2s + 1) which is not in general an integer, but 
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we can still write as = n + r] where n is an integer and < 77 < 1 . This eigenvalue is still 
consistent with Eq. (|1.2|) given in the Introduction for the area eigenvalues. Hence, we shall 
define a subspace Ti^ of Tif as follows: 

7Y^ = spanj I s„ = n + 77,7713)} (3.4) 

where rz = 0, 1, 2, ... and m^, = —Sn/^, —Snjl + 1, ..., s„/2. In the same manner, the third 
requirement on H implies that the time evolution of the BH takes place in ?i^. Hence, it 
is clear now that the third demand on the Hamiltonian is consistent with the spectrum of 
Eq. ( |1.2| ). From now on, for simplicity, we choose r/ = (note that 7i = 7irf=o). 

We have explained why the third demand on H is necessary in order to shift attention 
from the full Hilbert space to the subspace Ti which yields the degeneracy k^ for the number 
operator N. We shall now be interested in the kind of interaction term V which implements 
the third requirement. 

We introduce raising and lowering operators E^ and E of A, respectively, which we shall 
define explicitly later. These operators must satisfy the following commutation relation 

[A, E^] = E^ and [A, E] = -E. (3.5) 

From the third requirement on H it is quite clear that V contains E^ and E since (E)", for 
example, allows only transitions between states with hyperspins s and s — /c"/2. In addition, 
we are not interested in products of the form E'''E since they keep A unchanged (that is, 
[A, E'l'E] = 0). Combining all of this with the hermiticity of V we find that V can be written 
as 

V = h{A,E) + h*{A,E^) (3.6) 

where h is a complex function of A and E. 

Now we shall expand the function h a.s a. power series in E. That is, the first term 
is proportional to E, the second one to (E) etc. It is reasonable to assume that in a 
perturbative theory the transitions n ^ n±l have more weight than higher order transitions 
(such as n — > 72 ± 2). Hence we shall assume that up to a given approximation, h is linear 
in E. Later on, we shall justify this assumption by a comparison with the HR for a massive 
BH. But note that despite the linearization of h, the time-evolution operator, which is an 
exponential function of E and E^, contains all orders of these operators. We can thus write 
the "interaction" term as 

V = /(A)E + Etr(A), (3.7) 

where / is some function of A. If we normalize the raising and lowering operators to have 
expectation values of the order of unity (as we shall see shortly), it follows from the second 
requirement on the Hamiltonian that the function f{n) approaches zero as the parameter 
n {n = (A)) approaches infinity. The question is how fast f{n) approaches zero. We shall 
give here a motivation that f{n) ~ 1/v^ for large n, and in section ^ we shall estimate the 
function f{n) in a more deductive manner. 

It was shown by Kastrup [0 that there is a small correction to the BH entropy calcu- 
lated in canonical ensemble of the order of the logarithm of the area. Likewise, we [Q have 
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obtained a logarithmic correction to the entropy by using the grand-canonical approach (we 
have shown that the grand-canonical approach is also compatible with the HR). Padmanab- 
han [^], as well as Kim, Lee and Ji ||32|, have shown that for a small number of particles in 
the background of the Schwarzschild metric, (or if the system obeys Boltzmann statistics), 
the entropy is proportional to the logarithm of the area. This implies that there might be 
a correction to the entropy of the order of the logarithm of A. Motivated by the above 
arguments we shall write the entropy as: 

S=-^A + b\n{dA) (3.8) 

where b and d are constants. Thus we find a small correction to the energy U of the BH: 

b' 

dU = Tends = dM + --dM (3.9) 

where b' is some constant. Note that the correction to Hawking temperature Tbh, if there is 
one, is absorbed in b'. Hence, the correction to the energy U is proportional to the reciprocal 
of M. This is one reason for guessing that the function /(A) in Eq. ( ^.71) is proportional to 
1/vA for large (A) values. As we have mentioned earlier, we shall later derive this result 
in a more rigorous way. We now turn to define the operators E and E''" explicitly. 
The most general form of a raising operator a'^ in the Hilbert space 7i is given by 

a^ = E E E C-^,m> + ^,rn,){n,rm\ (3.10) 

n=0 mi=l m2=l 

where C^^^^ ^^^ complex numbers. The lowering operator a is defined by taking the Hermi- 
tian conjugate of Eq. (|3.10| ). Now, according to Eq. ( |2.22| ), in the observable Hilbert space 



Tiob, the raising operator a''" transforms to ha"''h. That is, the raising operator in Tiob is 
defined in terms of the observable states \n) instead of Eq. (|3.10|) . Hence, we are able to 
define E and E^ as the Susskind-Glogower operators p3| . 



E=Y,\n){n + l\ and E^ = ^ |n + l)(n|. (3.11) 

n=0 n=0 

It is clear from the definition that these operators satisfy relation ( |3.5| ). We can define also 
the standard creation and annihilation operators, a and a"!", respectively, as 

a=(N+l)i/2E^ at = Et(N + 1)1/2. (3^^2) 

It is clear that N = a^a, and 

oo 

[at,a] = ^|n)(n| = lob (3.13) 

ra=0 

where the lob indicates a unit operator in Tiob (notice that it is not a unit operator in Ti). To 
write the interaction term ( ^.7] ) with the operators E and E''' defined in Eq. ( |3.11D requires 
knowledge of the quantum phase problem. Hence, before we continue with the formalism, 
let us first review this problem in accordance with Lynch p3| . 
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C. The Quantum Phase Problem 

Classically, the electric field vector E{f,t), of an electromagnetic field contained in a 
cubic cavity of edge L, may be expanded in Fourier components: 

E{f, t)(xJ2 (flfce^^^-^^-^*) + a|e-^(^-^^-^*)) (3.14) 

k 

where we assumed for simplicity a linearly polarized field at a single frequency uj. The sum 
is taken over all the plane wave modes, where k is the wave vector. If we are limited to a 
single excited mode, say, k = ko, we can write 

E{r, t) oc (ae*('=^o-"^-^*) + a^e-^^'^^o"^"^*)) (3.15) 

where a^ a^ . Writing a = re^'^ we finally obtain 

E{f,t) = Eocos{ko-f-ujt + (j)). (3.16) 

In this way we decompose the classical electric field into amplitude and phase components. 
In quantum mechanics a ^ a and a* -^ a^ where a and a} satisfy the canonical commu- 
tation relation ( p.l3|) . Dirac (1927) was the first to attempt to define a phase operator. He 
decomposed a into an amplitude and a phase component 

a = e^^N^/2 (3.17) 

where N = a^a. The commutator [a, a''"] = 1 leads to the Lerner criterion 

[e'-^, N] = e^"^. (3.18) 



Expanding e**^ on both sides of relation (|3.18|) gives 



[N,0]=i (3.19) 

and thus shows that N and are canonically conjugate operators. This relation leads to 
the number — phase uncertainty relation of quantum optics: 

ANA0 > -. (3.20) 

That is, to create a radiation with a sharp phase, such as a laser beam, we would need to 
have many photons. Similarly, as we shall see later, a massive BH (analog of many quanta) 
has a sharp phase. 

But Eq. (|3.19|) leads to serious problems. For example, if we take its matrix elements in 
the number state basis, we find that 

{n — n')(n|(/)|n') = ibn^n' (3-21) 

which leads to the impossible equation =i for n = n' . Moreover, Susskind and Glo- 



gower |3^ have shown that e*'' is not unitary, or equivalently, (j) is not Hermitian. 
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The recent paper by Bojowald, Kastrup, Schramm and Strobl ||28| suggests a solution 



to these problems. However, as we shall see shortly, despite these problems, if one does not 
define the phase as an operator (just as time is not an operator in quantum mechanics), the 
above criticism is not a decisive objection to the phase- number uncertainty relation ( p. 201) . 



For although we have no time operator, we have time-energy uncertainty relation. The 
reason why it is difficult to define the phase or the time operator is the fact that the particle 
number or energy spectrum is non-negative and, therefore, is bounded from below. This 
situation is drastically different from the one for position and momentum whose eigenvalue 
spectra are unbounded. 

Let us now compare the Dirac definition for the phase operator Eq.( |3.17|) , and Eq. ( p. 121 ) 



One can see that E and E^ are analogous of e**^ and e"*"^ respectively. Now, using the 
definition Eq. ( p.ll| ) it is easy to show that E and E"'' are "one-sided unitary" or an isometry. 
That is, 

EEt = l but E^E = 1- |0)(0|. (3.22) 

Despite the fact that the vacuum projection "spoils" the unitarity of E, for states with a 
very large number occupations (that is, (N) ^ 1) and hence small vacuum component, we 
can treat E as an approximately unitary operator. As we shall see now, the last statement 
helps us to define the phase of a BH. 



D. The Phase of a Black Hole 

For a BH not near the Planck mass, (N) ^ 1. Thus we can define a phase operator by 

e"f' ^ E. (3.23) 



Hence, the interaction term V, which has the form of Eq. (|3.7| ), describes an interaction 
between the area operator A and the phase operator 0. As we shall see later, the HR in our 
model is due to an area-phase interaction. 

The definition of the quantum phase of a BH in Eq. ( |3.23|) raises the question, what is 
the physical interpretation of the "new" observable (j) ? Furthermore, what are the physical 
grounds for an area-phase interaction appearing in the Hamiltonian? In order to answer 
these questions, we first note that the phase operator is the canonical conjugate of the area 
operator (see Eq. p.l9| )). This is reminiscent of work by Carlip and Teitelboim 



Carlip and Teitelboim have shown by using the standard (Euclidean) action principle for 
gravitational field, that the horizon area and the opening angle at the horizon are canonical 
conjugates. They start with the wave functional %jj = ip[T,^^^ g] where the separation time 
at spatial infinity T, is an extra argument (in addition to the intrinsic geometry of a spatial 
section ^^^g) since space is not closed. The canonical conjugation relation of T and the ADM 
mass M is expressed in the Schrodinger equation 

h dip , , 

-^ + MV^ = (3.24) 

which can be obtained from an extended Wheeler-DeWitt equation (that is, the appropriate 
action is the canonical action Sm+g supplemented by —TM). In the same manner they 
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introduce an extra parameter ( "dimensionless internal time" ) associated with the horizon 
by the relation 

|||-.„^0 (3.5, 

where A is the horizon area. They work in a system of polar coordinates r,T in 3f?^ (the 
spacetime topology is 3fJ^ x S''^^^) where they set the origin r = r+ to be the horizon of the BH 
and the angle r to be the Killing time (for more details see Carlip and Teitelboim paper [^). 
Next, they show that the parameter iQ is " the total proper angle (proper length divided 
by proper radius) of an arc of very small radius and coordinate angular opening ri — T2" . 

The main result of Carlip and Teitelboim is that the canonical conjugate of the horizon 
area is the opening angle. The Euclidean continuation of the hyperbolic angle 6 is bounded 
and even periodic. This not only favors our simple model, but it gives a motivation for the 
assumption that the horizon area (in the Euclidean sector) can be represented by a number 
operator ! We are able to define the phase in Eq. (|3.23| ) since we have assumed that the area 
spectrum is equally spaced. If, for example, the area spacing a„+i — a„ approached zero for 
large quantum number n (as suggested by Ashtekar and others |^), the quantity canonically 
conjugate to area would not be a phase operator (its eigenvalues would not be bounded). 

In our simple model, the canonical conjugate of the horizon area is also some kind of an 
angle. Since in our model the phase (angle) is an operator, we shall identify the eigenvalues 
of with the "opening angles" at the horizon, as discussed by Carlip and Teitelboim. As 
we shall see below, a further analogy between our approach and theirs can be found in the 
description of the HR. 



Recently, Massar and Parentani [^ have shown that the probability for a BH to emit a 
particle is given by 

Pm^m-x = NiX, M)e-^^(^'^)/^ (3.26) 

where AA is the change in the horizon area induced by the emission, and A^ is some phase 
space factor. Since they did not neglect the specific heat of the BH, expression ( p.26|) 
improves Hawking's result. We shall now point out the origin of this result. 

As suggested by Carlip and Teitelboim, Massar and Parentani have supplemented the 
canonical action with two boundary terms, —TM (boundary at infinity) and QA/Sn (bound- 
ary at the horizon). Hence, they face four different actions according to which quantities are 
fixed in the variational principle. They fix the ADM mass M and the opening angle in 
order to describe a physical situation in which the BH and the surrounding matter exchange 
energy, but no energy is exchanged with infinity (this is also the process we are studying). 
It is through the boundary term Q A/Sir that they obtain Eq. (|3.26|) . 

In our approach, as we shall see later, the HR appears as a result of the interaction ( p.TD , 
an interaction between the phase ("opening angle") and the area of the BH. And Massar 
and Parentani have obtained Eq. ( p.26| ) after introducing the boundary term QA/Stt in the 



action; this can be considered (in the context of our model) as a 'coupling' between the 
horizon area and the opening angle. Hence the Massar-Parentani approach is compatible 
with our model. For this reason we identify the area-phase interaction term in Eq. ( p.7| ) 
with the boundary term QA/8tt. Note that in our approach the interaction term appears in 
the Hamiltonian, whereas in the Massar-Parentani approach the boundary term appears in 
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the action. Furthermore, it is interesting to note that the effect of general relativity (that 
is, curved spacetime) appears in our simple model as a coupling between the observable A 
and its own canonical conjugate ! 

We now ask what our model has to say about radiation by the BH. In addition to the 
Hamiltonian we need specify the state of the BH. 

IV. THE INITIAL STATE 

One might think that the initial state is just the an area eigenstate |no) where M = 
mQ^fn^ is the initial mass of the BH. But as we shall see in section [V|, an area eigenstate 
describes a static BH (that is, no HR). The question is then, what is the initial state which 
leads to the semi-classical radiation (HR), and why an area eigenstate is not an appropriate 
candidate. 

The question above is not unique to our problem, but is a general question in quantum 
mechanics: What are the states that most closely imitate classical behavior? The simplest 
case where this question appears is that of the simple harmonic oscillator. Hence we shall 
first review, very shortly, what are the states that imitates the classical harmonic oscillations 
(see Sakurai). 

In the Heisenberg picture the position operator x(t) of an harmonic oscillator is given by 

x(t) =x(0)costut+ I ^--^ ) sincjt (4.1) 

y muj J 

where uj is the angular frequency of the classical oscillator and p(0) is the momentum 
operator at t = 0. The expectation value (n|x(t)|n) vanishes because the operators x(0) 
and p(0) change n by ±1 and \n) and \n + 1) are orthogonal. This result is different from 
the classical result no matter how large n may be. In order to obtain the classical result, we 
have to take the expectation value with respect to a coherent state. 



Let us perform the Bogoliubov transformation for coherent states ||35|| : 



a^a(A) = a-A (4.2) 

where A is a c-number. Then, a coherent state is defined as the vacuum for q;(A), and is 
denoted by |0(A)): 

a(A)|0(A)) = or a|0(A)) = A|0(A)). (4.3) 

We then have 

a(A) = U,(A)aU;^(A) with Ue(A) = exp (zGc(A)) (4.4) 

in which the generator Gc(A) is given by: 

Ge(A) =i(A*a-aU). (4.5) 

We note parenthetically that if we allow A to be a function of A, we shall be able to express 
the interaction term in Eq. ( p.7|) as the generator of BH coherent states. That is, 
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V = Gc(A(A)) with X{A)=if*{A){N + l)-^/\ (4.6) 



This remarkable result implies that the interaction term given in Eq. (|3.71 ) induces the 
Bogoliubov transformation ( |4.2| ); we shall see that it defines the BH states which exhibit 
HR. 

Let us now list two properties of a coherent state which will be relevant to our problem. 
Firstly, a coherent state may be expressed as a superposition of the number (N) eigenstates 

|0(A)) = U,(A)|0) = 5:g(n)|n) (4.7) 

n=0 

where the distribution of |g'(?^)P with respect to n is of the Poisson type about some mean 
value n = |Ap: 

\qin)\'=^—jexp{-n). (4.8) 

Secondly, a coherent state satisfies the minimum uncertainty product relation. 

By analogy with a simple harmonic oscillator, we shall ask how can we construct a super- 
position of the BH area eigenstates that most closely imitates the semi-classical result (HR)? 
In comparison with the harmonic oscillator, would it be some kind of a coherent state? We 
shall call the state that does the desired job the coherent-termal stateQ and denote it by 
\CT; no) where no = {CT; no|N|CT; uq). Hence, in the limit of a massive BH we demand 

\im{CT;no\M(t)\CT;no) = (M(t)) (4.9) 

no>l 

where {M{t)) is the result obtained in the semi-classical approach (see section ^. 

An observer at infinity is limited by the energy-time (area-phase) uncertainty principle 
when he tries to determine the mass (area) of a BH formed by a gravitational collapse. If 
the BH is an eternal one, in which one can determine exactly its mass, there will be no HR 
because the initial state will be a mass eigenstate (see section |V|) . By contrast, if the BH 
has been formed by gravitational collapse, there is always radiation in the region between 
the BH and the observer. Even if hypothetically at some initial time, say to, there is no 
radiation, still the measurement takes time, say At, in which a small amount of radiation is 
produced according to the mass-time uncertainty principle. Hence, in order to be compatible 
with the HR, the initial state must satisfy ((AM)^) ^ 0, or equivalently, the HR is due to 
the uncertainty (of an observer at infinity) about the area (mass) of the BH. 

The coherent-thermal state (CT-state) may be expressed as follows: 

oo 

\CT;no) = Y.Pn^'^M-Hn))\n) (4.10) 

n=0 

where P„ is the distribution with respect to n and exp(— ^//(n)) is some phase coefficient 
where z/(n) is a real function of n. We can set exp(— «//(«)) = 1 by defining \n)^ew = 



^The termal feature of this state has been actually erased after performing the projection operation 
of the horizon. 
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exp{—ih'{n))\n) and doing all the previous analysis with respect to these new states. The 
new Susskind-Glogower's operators, E and E''', are different from the old ones, but the form 
of the Hamiltonian and the interaction term given in Eq. ( |3.7| ) are unchanged since /(A) 
can absorb the changes. Hence we conclude that the CT-state may be written as: 

oo 

\CT-n) = Y.Pn^»- (4.11) 

We are interested in the semi-classical limit uq -^ oo. In this limit (Jm/ (^) —>■ where 
cr|^ = (A'P) — (M)^. Hence, in the limit of large uq value, the distribution P„ is of the 
Gaussian type about the mean value no [^. That is, 

1 
/ a \2 / a , ,2 



where we have left the variance, a\ = 2no/a, to be determined later. But a distribution of 
the Poisson type (see Eq. (f4.8|) ) becomes, in the limit no — > oo, Gaussian with a variance 
cr^ = no. Thus, for a coherent state a = 2. However, to be more general we shall not yet 
set a = 2. 

Let us now compare the semi-classical results and our quantum mechanical approach 
when the averages are taken with respect to the CT-state 

\CT; no) = (-^) " £ exp (-^{n - no)^) \n). (4.13) 

V47rno/ ^Q V 8no / 



V. THE ROLE OF HAWKING RADIATION 

We start by emphasizing the similarity of the following two roads to the HR: 
(l)If one attributes to a BH an entropy Sbh = v4/4£p (Bekenstein [0), it leads to a 
temperature Tbh = h/SiiM (Hawking JT^); hence the BH radiates and decreases its horizon 



area {A = IGttM^). Due to the thermal character of the radiation, one can calculate the 
rate of change of the BH mass using Stefan-Boltzmann law (see Eq. (|5.1|) ). 
(2)BH entropy, together with the assumption that the horizon area is represented by the 
number operator A = 41nfc£pN, implies a huge degeneracy with respect to N, k^. This 
degeneracy motivated us to introduce into the Hamiltonian an interaction term of the form 
of Eq. (|3.6| ). A BH with such high degeneracy acts in order to reduce its internal degrees 
of freedom. Hence the area of the BH is not conserved, i.e. [H, A] 7^ 0. Later on we shall 
see that it is indeed possible to choose the interaction term such that the area of the BH 
decreases in time. 

The purpose of this section is to compare these two roads. Such a comparison examines 
the validity of our assumptions that lead to the form of the Hamiltonian given in the previous 
section. Moreover, we shall be able to derive the form of the function /(A) in Eq. (|3.7|) for 
large (A) values. We shall point out that these two roads should match only in the limit of 
large (A) values (the semi-classical limit). For much smaller values of (A) (near the Planck 
scale), our approach predicts new results. 



Due to Hawking radiance and its thermal character, it is possible to estimate the rate 
of change of the mass (M) of the BH. Using the Stefan-Boltzmann law for a surface area 
167r(M)^ at temperature Tbh = (^/8^)(^)~^5 we find (see Bekenstein and Mukhanov ^ 



and Page p5|] ) 

d{M) -ffi 



dt 153607r(M)2 (M)^ 



(5.1) 



where 7 lumps uncertainties about the species of particles emitted. We shall assume that 
the mass of the BH at time t = is equal to the total energy of the system, that is, (M(0)) 
is the ADM mass. Now, solving Eq. { p.l\) we find 

(MW> = (M(0)) l-(^^i^) =(M(0)>-^^^^<-(^^^<=-..- (5.2) 



where e is defined in Eq. (|5.1| ) and the mass of the BH is assumed to be very large compared 
to {etf. 

In our model, the mass of the BH is represented by the mass operator M = M^ = M/^(0) 
where the subscripts S and H stand for the Schrodinger and Heisenberg pictures, respec- 
tively. The evolution in time of the mass operator in the Heisenberg picture is determined 
by 

Mnit) = exp (^at) M(0) exp (-^Ht) = M(0) + ^[H, M(0)]t 

+ ^Q)'[H,[H,M(0)]]t^ + --- (5.3) 



where the Hamiltonian is given by Eq. (|3.3| ) and Eq. ( p.7|) . In order to compare Eq. ( ^.3| 
with the semi-classical result ( |5.2| ), we have to take the average of both sides of Eq. ( p. 3D 
We shall first compare between the coefficients that are proportional to t in both Eq. ( ^.2| 
and (the mean value of) Eq. (E 



A. The First Order Corrections 

The Hamiltonian ( p.3| ) may be written as 

00 
H = M(0) + J2 (/Hl^)(^ + 1| + riri)\n + l)(n|) (5.4) 

n=0 

where we used Eq. ( p.7|) for V. Hence it clear that 

[H,M(0)] = [V,M(0)] = 



mo J2 [iV^^+^ - v^) ifin)\n){n + 1| - r{n)\n + l){n\)) . (5.5) 

n=0 

Note that if we take the average with respect to an area eigenstate \n), we would find that 

(n|[H,M(0)]|n) = 0, (5.6) 
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and it would be impossible to compare with Eq. 



Now, taking the average with respect to the CT-state ( 4.13 ) we find 

oo 

([H,M(0)])cT;no = 2moE /(^)(V^^"TT- V^)(|n)(n+ l|)cT;no 

n=0 

2mo ^/H V^^^TI-v^eV «"o^ b„o^ -^ 



n=0 



'^"-^ 2mo/(no)(v/^^^TI- v^)e^^) (5.7) 

where we have assumed that fin) is purely imaginary since we want that (H)cT;no = 
(M(0))cT;no (the mass of the BH at time t = is equal to the total energy of the sys- 
tem). Now, comparing i^~^([H, M(0)])cT;no with the coefficient oft in Eq. ( |5.2| ) we find 

/K) "°=^"^ ,/f' , e(^) (5.8) 



where we have taken v^no + 1 — ^Jn^ ^ l/{2y/no) since no is very large. Because (M(0)) 
fno\/^ we conclude 



f^n)"^'%i^)^. (5.9) 

Note that the no in the exponent was not replaced by n because this exponent originally 
came from the initial state (see Eq. ( |5.7| ) before taking the limit of large tiq). This equation 
holds only for n ^ 1 and there are corrections in the limit of small n. Note that this result is 
in harmony with our previous estimate of /(A) as proportional to 1/vA (see the arguments 
below Eq. (gj)). 

It is clear that the function f{n) should not be dependent on the CT-state characterized 
by the number no (otherwise, the Hamiltonian would be dependent on the states we take 
average with). Hence, the factor exp [j^] should be independent of no- There are two 
ways to achieve this. One is to assume that a is proportional to uo such that the ratio a/no 
is independent of tiq. The other way is to assume that a is of the order of unity, and since 
our calculations are valid up to corrections of O(l/no), we can replace exp (j^) by one. 

By contrast, it would be incorrect to replace exp ij^) by unity if, for example, a ~ y/no- 
The reason is that we cannot neglect corrections of 0{l/y/no) and thus the function /(n) 
will be dependent on no by a small correction (~ 1/y/no) which cannot be neglected. Since 
the Hamiltonian should not depend on no (even not as a small correction), we conclude that 
either a is of the order of unity, or it is of the order of no. The question now is which one 
of the two possibilities for a is the correct one. 

One may ask if the area-phase uncertainty relation (|3.2CI|) can supply the required infor- 



mation about a. As we have mentioned earlier, a coherent state (e.g. the initial CT- state 
with a = 2) satisfies the minimum uncertainty product relation. In Appendix C, we show 
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that both of the choices for a lead to the ininiinuiii uncertainty product relation ( p. 201) . 
Thus we are left with an unknown parameter. 



As we shall see shortly, Eq. ( p.9| ) leads to a second order correction of the form 
~ t^/(M(0))^ for the evaporation of the black hole [see Eq. ( ^.2|) ]. This result is, sur- 
prisingly, consistent with the semi-classical result Eq. ( |5.2| ). Hence, we shall compare also 
the coefficients of second order in both approaches, and deduce restrictions on a. 

B. The Second Order Corrections 

In order to calculate the second order of (M(t)) with respect to the time t, we use 
Eq. ( |5.3|) and calculate the average of 

[H, [H, M(0)]] = [M(0), [V, M(0)]] + [V, [V, M(0)]] (5.10) 

with respect to the CT-state. Due to the identity (|n)(n + 1|) = (|t2 + 1)(?t-|), the first term 
in Eq. ( |5.1CI| ) averages out to zero. To calculate the average of the second term in Eq. ( ^.101) , 



we substitute Eq. ( [5. 81 ) in the expression for V. Thus, we find 

fc2 2 oo , / 1 T 

([V, [V,M(0)]])cT;no = -X«^^^ E -^ - 7— TWI 1 (I^)(^I)CT;.0 



m: 







1 




\n\/n + 1 (n + 1)V"- 



___e<=^>^3+E(|. + 2)(„|),,„.) 

t2 2 

'3e(s%)+e(-8^)^ (5.11) 



2(M(0)) 

where we have used the fact that the contribution of small n's (compared to rio) may be ne- 
glected under the average. Finally, comparing Eq. ( |5.11| ) with the semi-classical result ( p.2|) 
we find 

3e(8^) +e^"«^^ =4. (5.12) 



Note that Eq. ( [5. 121 ) has two solutions. One is exp(a/8r2o) = 1. This solution is plausible, up 



to a very good approximation, if a is of order of unity since then exp(a/8no) = 1 + O(l/no). 
A coherent state, that is a = 2, also satisfies Eq. (|5.12|) . This remarkable result shows 



that our model for the Stefan-Boltzmann law agrees with the semi-classical result up to 
the second order in time. We have used the first order correction to obtain /(n), and this 
allowed us to reproduce exactly the second order correction (including the non-dimensional 
numerical coefficient). 

The other solution of Eq. ( |5.12| ) is exp(a/8rao) = \ which is impossible for a > 0. Hence, 



it seems at ffist sight that a solution for a of order no is impossible. 

In the last argument we have assumed that the interaction term in Eq. ( |3.7| ) is exact. 
If one insists that the initial CT-state, which imitates the semi-classical result, is such that 
a ~ rio (and not of the order of unity), we must then add a correction to the Hamiltonian. 
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In appendix B we have estimated the second order correction (if there is one) to the Hamil- 
tonian, and we have reproduced Eq. ( ^.12 ) with a correction coming from the extra term 



in the Hamihonian. We shall emphasize here that the second order correction contributes 
Eq. ( |5.12|) because we have calculated (M(t)) up to a second order in time. 



VI. STATIC BLACK HOLES 



At first glance it seems impossible to describe a static BH in our model since [H, A] 7^ 0. 
But, as we have seen in section |VA , the first order correction to the rate of change of the 



mass of the BH is zero when the initial state is taken to be an area eigenstate. Thus, area 
eigenstates seem to be good candidates for describing eternal BHs. Note that these states 
satisfy ((AM)^) = which is compatible with eternal BHs since no gravitational collapse is 
involved (see section fV\) . As we shall see in this section, for a BH (not a primordial one) 
which is described by an area eigenstate, we have to wait a long time (much greater than 
the age of the universe) in order to observe any change in its horizon area, and even then it 
will only be a small fiuctuation about the mean area value. 

The idea that area eigenstates describe static BHs has interesting implications. For 
example, if one had a 'device' able to measure the area of the BH (that is, the BH state 
collapse into an area eigenstate), it would, by measuring it, stop the time evolution of the 
BH! This means that our knowledge about the area of the BH plays a central role, and 
affects its mass evaporation. 

We start with the Hamiltonian given in Eq. (^.4|) where f{n) is given by Eq. ( |5.9| ) for 
large n. As we have done in the previous sections, since f{n) approaches zero for large n 
values, we shall treat the interaction term in a perturbative manner. This time we shall 
work in the interaction picture. 

The time dependent "potential" in the interaction picture is given by 

Yj{t) = exp (^M{0)t] Vexp (-^M{0)t] (6.1) 

where V is the interaction term given in Eq. ( ^.41) . The time evolution operator (in the 
interaction picture) is then expressed as 



U/(t) = T exp (-th-^ I Yi{t')di 

= 1 _ 1 /■* Y,{t')dt' + {-'-f f dt' C dt"Yi{t')\i{t") 
n Jo n Jo Jo 



(6.2) 



where T is the time ordering operator. We proceed to evaluate the average of the area 
operator when the initial state is the area eigenstate |no). We shall calculate the average up 
to the first non vanishing order. That is, 

{A)t = {no\V\{t)AUi{t)\no) = n^ao 

+ ^2(^0! fYi{t')dt' A fdt"Yi{t")\n,) 
n Jo Jo 

- ^(nol Jl dt' f^ dt"Vi{t')Vi{t")\no) + O (V,^) (6.3) 
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where we have used Eq. ( |6.2| ) (note that the area operator in the interaction picture is the 
same as in Schrodinger picture). 

Now, using Eq. ( |6.1| ) and Eq. ( [5.4| ) to express V/(t), we find after some algebraic ma- 
nipulations that Eq. ( |6.3|) can be written as 

(A)i = nooo 

4ao (\f{no-l)l . , fl \ |/K)p . , (I \\ 



(6.4) 



where 



Wno,no±l = T = T ■ \0-^) 

Now, because f{no), and cUng^^o-i-i are of O(l/iyno) for uq ^ 1, we find that |/(no — 



l)P/a;^^ „g_]^ is equal to |/(^o)P/^.^o "o+i ^P ^° ^ ^^^^ small correction of the order of I/uq. 
Thus, after using Eq. (|5.9|) for f{n) we find 

for Kq ^ 1. Substituting Eq. ( |6.6| ) in Eq. ( |6.4| ) we finally obtain 

{A)t = aoTlo - ao^^ fsin^ (-UJno,no^lt] - Sin^ (-UJno,no+ltj] 

= aoUo — ao^^ sin(co'''"'^^t) sin(ti;*-^-'t) (6.7) 

where (with G and c displayed) 

,(1) - _!!ii /'£!'! , ,(2) - ('"*)'* r <>' ^ .... . _ "-"^"f (s;) 



-" = 8^ Ig j • '^" = 118^ b^j -'' « = ^400^- («^«' 



Note that ^ is a dimensionless coefficient, whereas 7 was first introduced in Eq. ( |5.1| ). 

The frequency cu^^-* does not contain h; hence we call it the classical frequency of the BH. 
The second frequency, uj^'^\ is the quantum correction to uj^^\ and even for a microscopic BH 
with (M(0)) = 10-^^Kg, co'(^) ^ 10~^^sec~^ (where we assumed InA; ~ 1). Thus sin(co'(^)t) ?» 
for t not much larger that the age of the universe. We conclude that a BH (not a primordial 
one) which is described initially by a pure state \no) will stay in this pure state since {A)t = 
qqu over the lifetime of the universe. 

The frequency cu*^^^ becomes significant when l/cj*^^^ is smaller then the age of the universe. 
That is, BHs with a mass smaller than lO^^Kg. For these BHs it is impossible to neglect uj'^'^\ 
But still, if we assume that ^ (or 7) is not extremely large, Eq. ( |6.7|) describes a basically 
static BH with small time dependent fiuctuations around its mean area value. 
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VII. SUMMARY AND CONCLUSIONS 

We have described in this paper a fully quantum mechanical model for a Schwarzschild 
BH. Starting from the assumption that the horizon area is quantized with equally spaced 
spectrum, we have obtained for the black hole observables an algebra of angular momentum 
(the hyperspin of a Schwarzschild BH) with the restriction that the hyperspin can only have 
the values s„ = (A;" — 1)/2. This serves as motivation for the introduction in the Hamiltonian 
of an interaction term. Surprisingly, this term leads to the correct properties of the HR. 

We have introduced the phase of a BH as the canonical conjugate to the area operator, 
and identified it with the "opening angle" of Carlip and Teitelboim E^ . Then the "interac- 



tion" term in the Hamiltonian has been interpreted as a coupling between the horizon area 
and the phase of the BH. This coupling we associated with the boundary term QA/Stt in 
the action of a BH spacetime, which according to Massar and Parentani |3^ corresponds to 
the HR. 

To discuss the HR it was required to choose an appropriate initial state which most closely 
imitates the semi-classical result (HR). We have chosen the initial state to be Gaussian 
distributed with area about some mean value hq. Hence, the variance a^ of this state was 
left undetermined. However, after the comparison with the HR, it has been shown that 
either ^^/(A) ~ Mp/M, or aA/{A) ~ {Mp/Mf. 

For the first possibility, the Hamiltonian which describes the system is exactly of the 
form of Eq. ( |5.4| ). On the other hand, if aA_/{A.) ~ (Adp/M)"^, then the Hamiltonian given 
by Eq. (|5.4| ) is just part (in a perturbation theory) of the full Hamiltonian of the BH (see 
Appendix B). 

We have seen in the previous section that an area eigenstate describes a static BH. Hence, 
in our model the HR is due to the dispersion in the area eigenstates. It was shown long 



ago by Bekenstein [^ that the fluctuations in the mass of the BH, and thus also in its 



area, increases with time. Also Wu and Ford [^ have shown recently that the variance of 



the mass grows linearly in time. It is then clear that as long as the BH looses mass, the 
observers loose their knowledge about the area of the BH. Perhaps this is the way that a 
BH avoids a singularity at the culmination of its evaporation. 
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APPENDIX A: SETTING [/(N) = 1 



In order to prove that n(N) = 1 in Eq. (|2.6| ), we first define three operators in the same 
manner as in Eq. (2.14), that is 



Ji^i(/(N)g + gtr(N)) 
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j. = ^(/(N)g-gtr(N)) 

J3 = |/(N)r(G-V-i) (Al) 

where m(N) = |/(N)p since m(N) must have real positive eigenvalues (see Eq. ( p.6|) ). These 
three operators satisfy the angular momentum commutation rules and thus the eigenvalues 
of J3 and J^ are m^ = — j, — j + 1, ..., j and j{j + 1), respectively, where j is an integer or 
half integer. 

Using Eq. ([Al|) and Eq. ( p.6|) one can obtain 

4J2 = |/(N)|2 fe^ - 1 - (1 - |/(N)|2) (2G -k^ - if) . (A2) 

Hence it is clear that the eigenstates \n,m) of N and G are also eigenstates of J3 and J^. 
Thus we shall write m^ and j in terms of n and m as follows: 

2m3 = \l{n)\^ {2m - k" - 1) (A3) 

4j(j + 1) = |/(n)p (F" - 1 - (1 - \l{n)\') (2m - fc" - 1)^) . (A4) 

From Eq. ( [A3D it follows that |/(n)p is an integer since 2m^ is an integer. Thus, |/(n)p > 1. 
Setting for example m = 1 in Eq. ( [A3| ) and Eq. ( |A^ ) we find that 

2m3 = |/(n)p(l-A;") (A5) 

4j(j + 1) = |/(n)P (P - 1) (P + 2 - \l{n)\') . (A6) 

Now, the eigenvalues of an angular momentum satisfy the condition 2j > |2r7i3| which leads 
to 4j(j + 1) > |2r7i3| (|2r7i3| + 2). Hence, using both Eq. ([A 51) and Eq. (|A6|) we finally obtain 



|/(n)p (fc" - 1) (A;" + 2 - |/(n)|2) > |/(n)|2 (A;" - 1) (A;" + 1) (A7) 

which is satisfied only if |/(^)P < 1. Thus we have proved that u{n) = \l{n)\'^ = 1. 

APPENDIX B: A SECOND ORDER CORRECTION TO THE HAMILTONIAN 

The second order correction is related to transitions n ^ n ±2 (see section [111B| ), and 



hence may be written as: 

00 
H(2) = Y, ir{n)\n + 2){n\+r*{n)\n){n + 2\) (Bl) 



n=0 



where r(n) is a complex function of n. Now, since to zeroth order H*^"^ = M ~ vA, and 
the first order correction H*^^) = V ~ 1/vA for large (A) values, we conclude that for large 
(A) values the second order correction is proportional to 1/A^/^, i.e. each term of a higher 
order in the perturbation theory is smaller by one extra power of A. Thus, the function 
r{n) in Eq. (pl| ) can be written for large n as 
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■M'j^) 1 



rH = C^e^-o^-^ (B2) 



rrin n^/^ 



where C is a dimensionless number. It is reasonable to assume that ( is of the order of unity 
since the dimensionless coefficient of r{n)/mo in a perturbative theory is supposed to be of 
the same order as the square of the coefficient of f{n)/mQ. We might also assume that C, is 
purely real (and also negative), but this is not necessary because the imaginary part of C, (if 
it exists) will not contribute to Eq. (|5.12|) . 



We shall now take into account H*-^-* in the expression [H, [H, M(0)]] in Eq. ( |5.1CI| ). The 
only term in ([H, [H, M(0)]]) which contains H^^-', and is of the order of 1/(M(0))^, is 
([M(0), [H(2),M(0)]]). All the other terms are smaller by a factor of 1/(M(0))2, and may 
thus be neglected. A straightforward calculation gives 

([M(0), [H(2),M(0)]])cT;no = 2C^^[^e-^. (B3) 

Hence, we are now able to rewrite Eq. (|5.12|) , including the second order correction of the 
Hamiltonian, as 

3e(8%) + (1 + 4C)e^"^^ = 4 (B4) 



where C, is still to be determined. Note that if we choose C = (that is, H*^^^ = 0) Eq. ( P4|) 
reduces to Eq. ( p.l2| ), as it was expected. Thus, we conclude that if a is of the order of 



unity, then Eq. ( 5.12|) is satisfied automatically and we do not need a correction. This means 



that the Hamiltonian given in Eq. ( |3.3| ) with the interaction term ( p. 7] ) is equivalent to the 
exact Hamiltonian, at least up to O(/i^/(M(0))^). On the other hand, if a is 0(no), then 
the interaction term (|3.7| ) is only the first order term in a given perturbation theory. 

In order to determine a in this last case, we use Eq. ( p4|) . There are two solutions of 
Eq. (P4|) , but only one of them is plausible: 



i^) 2 + VI - 12/ 



&^' = . (B5) 



Note that / < since a > as it was expected. Relation (p5|) does not determine a, but 
instead express it as a function of the parameter /. 

APPENDIX C: SATURATION OF THE AREA-PHASE UNCERTAINTY 

RELATION 



We shall prove in this section that the CT-states defined in Eq. ( p.20|) saturate the 
area-phase uncertainty relation: 

ANA</>=^. (CI) 

The variance of N for the Gaussian CT-states is given by 



AN = J^. (C2) 

V a 
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We now turn to calculate A0. 

First we define sin in terms of the Susskind-Glogower's operators: 

sin0 = ^(E-ET). (C3) 

It is simple to verify that (sin (p)cT;no = 0; it is thus consistent to assume that also {4>)cT;no = 
0, where (p has been chosen to have its eigenvalues in the interval [— 7r,7r]. If a ~ 1, the 
variance AN ~ y/no is quite large, which implies that A0 is quite small. That is, is highly 
peaked around = 0. Hence we can assume that ~ sin0 which leads by Eq. ( |4.13| ) to 

i^'}cr:.. « ii.n.^f}cr:„ = i (l - exp (-£-)) « ^ (C4) 

for a ~ 1. Hence, 

and thus Eq. (|C1|) is satisfied. 



For a ~ no it is impossible to assume that is highly peaked around zero; instead we 
shall assume that the main contribution to (0^)cT;no comes from the region where —it/2 < 
< 7r/2. In this interval can be expressed as: 



(,2= lini '^^ 



where 

m 18 

= Y.Ci sin^' = sin^ + - sin^ + — sin^ + ... + C^ sin^™ 0. (C7) 

1=1 3 45 

The CT-states has the feature that 

m)cT;n, = {m)cT;n, = eXp (^-^) (C8) 

where / is an integer and a* = a/riQ is of the order of unity. Hence, expressing (0^)^ as 
polynomials in E and E^ and then taking the average according to Eq. ( P8|) , we were able 
to show (numerically) that 

((0^)jcT;no = f + 0((aT^^) (C9) 

were we have expressed the exponential functions in Eq. ( |C8|) as a Taylor series of a* . Hence, 
we conclude that for a* < 1 

which leads to Eq. (|Cl|) . 
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